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Abstract 


The entire field of astronomy is now in very rapid flux, and at 
the center of interest are problems relating to the very dense, rotat- 
ing, neutron stars now observed as pulsars, the hypothesized collapsed 
remains of stars known as black holes, and the not yet understood 
quasars. An understanding of these very exotic objects will undoubt- 
edly require the use of rapidly developing observational tools such 
as infrared astronomy, ultraviolet astronomy, and x-ray astronomy, as 
well as the classical tools of optical and radio astronomy. On the 
theoretical side it will be necessary to understand the behavior of 
very intense gravitational fields, such as are expected to occur in 
these objects, and the behavior of matter and radiation in such 
fields. 

We have used the so-called degenerate metric form, or Kerr-Schild 
metric form, to study several problems related to intense gravitational 
fields. These are as follows: 

I. The fundamental problem of relativistic gravitational theory 
is to solve the Einstein equations to obtain metrics which generalize 
the Lorentz metric of special relativity, 





(special relativity) 


We have worked with "algebraically special" space-times in order to 
simplify the mathematical problem. The classification of space-times 
is therefore treated first, and a new relation between the Petrov and 
Debever-Penrose schemes is presented. 


II. Our approach to the mathematical problem is to consider 
metrics of the specialized form 

v — Z/m. JL t JL , 

where £ is a null 4-vector, £^2. 11 = 0, and m is an arbitrary param- 
eter corresponding to the mass of the source of the field. We have 

studied time dependent metrics of this form and obtained very simpli- 
fied equations and a number of solutions, including the Kerr metric 
which describes a spinning black hole. Although the algebraic problem 
is solved in nearly complete generality, several special new cases 
arise which appear to describe gravitational waves, and are being 
studied further. 

III. The black hole solutions which arise in I. describe the 

exterior field of the final asymptotic state of a spinning star at 
the end of its evolution, if the mass exceeds a critical mass of 

about 2 solar masses. If the mass of a star is less than the crit- 

ical mass, it may instead become a pulsar — a rapidly rotating neu- 
tron star of about the same density as an atomic nucleus, -10 14 gm/cm 3 . 

Such neutron stars are incredibly complicated in structure and 

extremely interesting, especially those with a mass near the critical 
mass. These stars have an exterior field that is given, to first 

order in the stellar rotation rate, by the solutions discussed in I. 

We have studied the interior field for an idealized model, which is 
tractable analytically yet possesses some of the most interesting 
properties of very complicated neutron star models. The model is 
composed of an incompressible fluid undergoing a small amount of 


-ii- 



differential rotation. It is attractive in its simplicity and has 
properties such as a constant spatial density which are very close 
to much more cumbersome models. 
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1 . Introduction . 

In the decade following the publication of the field equations of 
general relativity^, attention focussed on methods of solution 
that exploited one or other of two simplifying assumptions: sym- 

metries in space-time, and linearization. This led to the 
Schwarzschild solutions ^ , the Weyl solutions the Einstein, 

De Sitter, Friedmann and other cosmological solutions and model 
universes ^ , and to first-order corrections to Newtonian gravi- 
tational theory^. Much of the emphasis was on the mechanical 
implications of the theory, in the sense that generalizations were 
sought of Newtonian mechanics, including the conservation laws for 
momentum, angular momentum and energy. 

In retrospect, an emphasis on analogies to classical mechanics 
looks to be too narrow. This point was emphasized in the late 

fc ?] 

1950’s by Trautman, Pirani and others v * . Instead of looking 

for a mechanical analogy, these authors argued that we should con- 
centrate attention on analogies with other systems of field equa- 
tions. In particular the electromagnetic field, which we know to 

(o') 

be intimately related to the gravitational field v , should be 
studied in a general Riemannian space-time. 

This approach was introduced by Pirani, Trautman and Lichnerowicz ^ , and 



2 . 


developed further by Debever^^ , Penrose^ Robinson^*^ and 
who were led to study the eigenbi vectors of the Riemann tensor as the 
natural analogue of the eigenvectors of the energy-momentum tensor of the 
electromagnetic field. Using the Einstein field equations, the Riemann 
tensor was found to define four null vectors at a general point of space- 
time These four null vectors, termed the Debever-Penrose princ- 

ipal null vectors, can be related to a coordinate invariant algebraic 
classification of space-tines developed by Petrov At any point of 

space- time, the space is defined to be algebraically /general if the four 
Debever-Penrose null vectors are independent, and algebraically special 
if two or more of the Debever-Penrose null vectors coincide. We can regard 
the study of algebraically special spaces as being equivalent to a simplify- 
ing assumption in the search for solutions of the field equations, just as 
physical symmetries and linearization are also simplifying assumptions. 


The first studies of algebraically special space-times emphasized applic- 
ations to gravitational optics and to gravitational radiation theory^ ^>16,17) 
However, the whole subject received a major new impetus in 1963, when Kerr'' y 
discovered a new solution of the field equations of great physical interest 
and importance. The Kerr solution, which appears to represent the exterior 
metric for the end-point of gravitational collapse of a massive rotating 
body (a "rotating black hole"^*^) has an algebraically special space-time 
structure and was found via an investigation of such space- times. Thus 
the latter were regarded with much greater interest and a search for other 
physically meaningful solutions of algebraically special type became logic- 
al. In the late 1960's, interest in the problem was again increased by the 
discovery of observational evidence that could be consistently interpreted 



(20 21 ) 

as arising from rota tine black holes ' . The evidence still admits 

this interpretation, and additional supporting observations were made in 


1972 


( 22 ) 


The main mathematical technique used so far to explore the structure of 
algebraically special space-times is differential geometry, couched in the 
tetrad formalism^ 23 ’ 24 ^ or in the spinor formalism ^ 25 ’ 26 ^ . These have the 
advantages of generality and power, and the disadvantage that they are still 
outside the mathematical repertoire of many physicists and astronomers. 

The plan in this work will be to study algebraically special space-times, 
using wherever possible and reasonable the most elementary and widely 
intelligible formulation. To accomplish this, we will find it necessary 
to compare and relate several alternative approaches, and in some cases to 
provide a completely new and elementary treatment. In Chapter 3 we disc- 
uss the classification of space- times, and in Chapters 4 and 5 we develop 
the relationship between two different methods of classification. After 
a short general discussion in Chapter 6, in Chapters 7 through 13 we pres- 
ent an elementary approach and solutions for a wide class bf algebraically 
special space- times. In Chapter 14» we briefly discuss rotating solutions 
in the presence of matter. 

Much of the material given here has been recently published or will 

soon be published in the Journal of Mathematical Physics, the Physical 

Review, and the Astrophysical Journal. Chapter 3, most of Chapter 4, 

and Chapter 6 are restatements of well-known material. The remaining 

work, except where references are quoted explicitly, are believed to 

(07 28 29 30 3ll 

represent new developments by the authors and co-authors v * ’ ’ * J . 
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2. Conventions and notation. 


We work throughout in a space-time of signature and with units 

such that c = G = 1 . 

In general, lower case Homan indices run from 1 to 3» upper case Roman 
indices run from 1 to 2 when associated with spinor quantities, and from 
1 to 6 when associated with matrices; and Greek indices run from 0 to 3. 
Unless otherwise explicitly indicated, the summation convention applies to 
all indices. 


We form the Ricci tensor by contraction on the first and third indices of 
the Riemann tensor, R^ v - • Covariant derivatives Eire denoted by 

a double line, k^^ , and ordinary derivatives by a single line, k^ly . 

Square brackets denote antisymmetrization, A \ “ A^, ) 

and round brackets denote symmetrization, A + A^w )• 

We associate a factor of l/N! with symmetrized or antisymmetrized sums over 
N indices. 


In general, the notation followed is that employed in Reference 8. 
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3. The Petrov Classification of Space-Times . 

In 1954, Petrov published a method of classifying space-times in terms of 
algebraic properties at a single point^^. This classification, which we 
will later show to be coordinate-independent, was made on a purely mathematic- 
al basis but was later found to be closely related to the physical problem of 
gravitational radiation (l5#17)^ 


The classification is based on the properties of the Riemann tensor in the 
vacuum case, and on the associated Weyl tensor in the non-vacuum case. For 
simplicity, we carry out the analysis for the case of a zero matter tensor, 
merely remarking that the non-vacuum case is handled in an exactly similar 
manner using the Weyl tensor in place of the Riemann tensor. 


Consider the Riemann tensor » which is antisymmetric inyu»y> and 

in , and is symmetric in the pairs • We ma P the 

indices onto the indices A, B by the homomorphism 

— ' 

yKV : 2.3 3 1 \% lo 2.o 3o / j 

A : l X 3 h S' ^ J 


Considered in terms of A and 5, Rr«r Is a second order symmetric tensor 

rt B 

with indices running from 1 to 6, and it can therefore be written as a 6x6 
matrix of the form: 



(3.2) 


where M and Q are real 3x3 symmetric matrices and N is a real 3x3 traceless 
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matrix (we assume the use of real coordinates). The form of (3.2) follows 
directly from the symmetries of the Riemann tensor and is isomorphic 

to the tensor • 

Raising the first index onT-A« and employing a local Lorentz coordinate 


frame leads to the form: 


-p A [A N/ 

K 5 ~ 

-^ T -q 


(3.3) 


for a zero matter tensor, the field equations K-v = 0 lead to the 


additional relations 


Q = - n 
N = N T 
'U(m) = 0 


(3.4) 


In this case we therefore have the form 


p n 


M N 

-ft M 


(3.5) 


where M and N are both symmetric real traceless 3x3 matrices. We now 

( 8 ) 

introduce the dual of the Riemann tensor, defined by v ' : 

* (3 - 6) 

where £ ^ is the usual unit tensor density^®^. We observe that mapping 

onto A, B as before and applying the field- equations yields the 


form for the dual tensor 


* f A _ 

N I - 


Is/ -tA 

n n 


(3.7) 



7. 


This result is most readily derived by element-by-element evaluation in 
the local Loren tz frame. We note here that the symmetries of the Riemann 
tensor lead to *R^ = 0 ’ whereas the symmetries of (which are the 

same as those of ) follow from = 0 . 

Thus, now defining the self-dual tensor K below, we have: 


A _ -pfl . • *,pA 




S 




n+wN -w(m+vN) 
Wm+ca/) 

This can be conveniently written in the direct product form: 


(3.8) 


-D+ ? -i? 

' i ? 




= ? ® r 


where P = M + i N 
and where 


J = 


“ I -c n 


(3.9) 

( 3 . 10 ) 
(3.11) 


P is a complex 3x3 traceless symmetric matrix and has 10 algebraically 
independent components. 

The eigenvectors of JC J are the direct products of the eigenvectors of 
P and J, and the eigenvalues . are the algebraic products of those of P and J. 
For J, we have eigenvalues 0 and 2 , and the corresponding eigenvectors 
(l, i) and (l, -i). Thus has at least three zero eigenvalues. For 

the Petrov classification, it is necessary to consider only the eigenvalues 
and eigenvectors of that are constructed using the non-zero eigenval- 
ue of J. These eigenvalues of may or may not be zero. Further, 

since the Petrov classification is in terms of invariant properties of the 
eigenvalues and eigenvectors of , the analysis applies in general 
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even though the particular form given by (3.8) holds only in the special 
(Lorentz) frame. 

The Petrov classification is performed according to the following simple set 
of statements about the eigenvalues and eigenvectors of the matrix P :- 

(a) If P has three distinct eigenvectors and three distinct eigenvalues, 
then space-time is Petrov Type I (this clearly applies only to the 
particular point considered - the Petrov Type may vary from one part 
of space-time to another). 

(b) If P has three distinct eigenvectors and two eigenvalues equal, then 
space-time is Petrov Type I-D . 

(c) If P has just two distinct eigenvectors and two distinct eigen- 
values, then space-time is Petrov Type II . 

(d) If P has just two distinct eigenvectors and has two equal (and 
hence zero, since P is traceless) eigenvalues, then space-time 
is Petrov Type II-N . 

(e) If P has only one eigenvector, then space- time is Petrov Type III. 

These properties of P are most readily established by reducing P to Jordan 
( 32 ) 

canonical form , and this suggests that we should also state the Petrov 
classification in an equivalent way in terms of the elementary divisors of 
P, as follows:- 

(a) If P has linear elementary divisors and three distinct eigenvalues 
then space- time is Petrov Type I. 

(b) If P has linear elementary divisors and just two distinct eigenvalues 



y. 


then space-time is Petrov TVT? 0 

(c) If P has just one linear elementary divisor and has two distinct 
eigenvalues then space-time is Petrov Type II. 

(d) If P has just one elementary divisor and equal eigenvalues 
then space-time is Petrov Type II-N. 

(e) If P has no linear elementary divisors then space-time is Petrov 
Type III. 


Finally, the Petrov classification is sometimes stated in terms of the 

(32) 

Segre characteristic of the matrix P, as follows 


Type: I 

I- D 

II 

II- N 

III 


Segre characteristic : ( 1 , 1 , 1 ) 

(( 1 , 1 ), 1 ) 
( 2 , 1 ) 
(( 2 , 1 )) 
(3) 


As we mentioned earlier, the classification is independent of the coordin- 
ate system. For, consider the coordinate transformation: x — ► 5c, so 



^ p/ 3 ' X 

30 


( 3 . 12 ) 



the coordinate transformation of (3.12) can be written as a matrix trans- 
formation of R 






(3.13) 
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where T*. = ^ — and T°, 

P*/ Px'' 

We can write this as the matrix product 


35c* ^ 


(3.14) 


H = T R T (3.15) 

But Tt = T\ f l ( ; ^ (5.16) 

Thus T = T ^ and so T is a similarity transformation on R. However, the 
Jordan canonical form, the eigenvalues and the number of linear elementary 
divisors are all unchanged under a similarity transformation, thus the 
Petrov classification of the space-time is invariant under coordinate trans- 
formations. 


The 3x3 complex matrix P is the fundamental matrix of the Petrov classification. 
We show in the next chapter that P is related by a similarity transformation 
to another 3x3 complex matrix arising from another (spinor) approach, and 
thus either matrix can be used as the basis for the classification. 

Although any coordinate independent representation of aspects of space-time is 
of interest, to this point the Petrov classification has appeared as a purely 
algebraic formalism, unrelated to the basic geometrical features of a space- 
time. In the next chapter, we introduce an alternative formulation in terms 


of spinors which leads to a direct and intimate connection with the space- 
time geometry. 



11 . 


4. The spinor approach and its relation to the Petrov classif ications . 

The Debever-Penrose principal null directions can also be used to classify 
space-times and in Chapter 5 we will do this explicitly. In this chapter, 
we wish to establish a spinor formalism and to relate certain fundamental 
spinor quantities to the Petrov classification matrix P. The Debever-Penrose 
principal null directions and their relationship to the Petrov classification 
can also be established by other methods^ , but the spinor approach 
appears to be the most revealing and the most elegant. 

The notation we will use is that of References 11, 26 and 34, and we will also 
quote and use a number of theorems concerning spinors without offering proofs, 
which can all be found in the cited references. 


Tensor and spinor quantities are transformed to each other via the general 
relationships: 


- | " 6 £ _ . — — ^ 

ci ~ I y 

— I ^ __ , -r-rti >— ci 

I v? ~ & AS I ci 

where the quantities satisfy the equations: 

fW + ftW* - fa s ** 


(4.1) 


(4.2) 


The S relate spinor space to tensor space, and Z. is a skew- 
symmetric metric spinor for the two-dimensional complex spinor space, enabl- 
ing spinor indices to be raised and lowered. Spinor indices take on the values 
1 and 2. Equations (4.1 ) extend in the obvious way to any number of tensor 
indices. In a tangent flat space with signature (+,-,-,-), a suitable set of 
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O’ ' s consists of multiples of the Pauli spin matrices, plus a mult- 
iple of the 2x2 identity matrix. 

To any second rank tensor F„ v there corresponds a 4-index spinor F A 
If i s skew, P^cj) can *> e she™ to have the unique spinor decomp- 

osition: 

^ABci> ^ £&j> + (4.3) 

where and fix are symmetric in their spinor indices. Further, 

if P^ v is real it may be shown that (Note that (4.3) 

differs from the notation of Reference 11 by a factor of two). 

If we apply (4.3) to the Riemann tensor R^-^v , which is skew in both 

r r A 

pairs of indices yUV and » the corresponding 8-index spinor may be 

shown to have the form: 

1 W* ^ ^AEPPCtf-DH ” (^"AScjt £<rH 

+ £c® £ef + £ftB fifcX £g-H + ^A-s£tj (4.4) 

The symmetry of in the pairs and p'X leads to the symmetries 

on X : 


X 


s %■ 


■ABc 2> 5 ^BAt* * 

In the vacuum case, applying the field equations leads to the additional 
conditions: 

^ ASci) ~~ ^ y Ascx C ^Xpnr fi (4.6) 


* 71. 


(4.5) 




Thus X is totally symmetric in all its indices. This important property 

( 1 ?) 

was apparently first remarked by Robinson . However, it was Penrose who 



first exploited the full significance of this fact, in making it basic to 


the classification of space-times, in a manner that we will shortly 
display^ As Penrose remarks, it is curious and remarkable that the 
complete symmetry of 9C follows only in a space with the signature 

Using (4.6), equation (4.4) reduces to the form: 

^ Xa£BFC£3>H ^ 

( y~Atc J> Z £f £c4 , + ZfU £ ca %££££ ^ 

The dual tensor *R can similarly be shown to have the spinor 
equivalent form: 

< > 


* ^ ^ V * 


(4.8) 


" * l ( Zepibu £ab ) 

Thus the self -dual tensor R^.v^o'X has the spinor equivalent form: 

= <2'X.M< 3 fep f(fii <4-9) 

Using (4.l) we then have explicitly: 

V” U* ASe , £er f«) (4 -’ o) 

or, raising the and v indices, 

*% - «- " (5 £«£«) ...... 

Since R + /* v ^ is antisymmetric in and in ^ we write (4.1l) in 
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the explicitly antisymmetrized form: 


^ A T^asch 


r* 


y CJ> 2 Jr v A 2 y ' 

A^a A c$ l,^ 

where in (4. 12) we define and ^ 

z 


cj> 


(4.12) 


r“> CD 




to be 


L 


/ 4V 

A 3 

- fe 

(^v/ - 

0 A 6 


■ fe 


c£ 


( ~ 

67, ^ 


(4.13) 


m\ - 
) £<£« 


Now we again introduce the tensor index mapping “ft , JT 

as in equation (3.1 )• Bars distinguish matrix and spinor indices and we have 

from (4.12) 


R**f - l 


1C 


AB 


(4.14) 


Y c.j> 

AB " c * A B 

Now, from equation (3.8) we note that the upper left 3x3 submatrix of £ £ 

is the matrix P. Thus if we restrict the range of the indices A and 
to 1,2 and 3, equation (4.I4) gives: 


»p A 

' 5 


yZ ^ab y c * 
A AS ^ *-2> A 5 


AS /vr c-2> 3 

where in (4.15) A and B take on the values 1,2 and 3, and A,B,C and U 
take on the values 1 and 2. 


(4.15) 


Equation (4.15) is an equation relating a 3x3 array, P, to a 4x4 array, X- • 
In general, the 4 - index spinor will have four eigenspinors, whereas the 
matrix P can have at most three eigenvectors. Thus it would seem at first 
sight that is not directly relatable to the Petrov classif iaction given 
earlier and based on the properties of P. However, since X symmetric in 
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all its spinor indices, any antisymmetric 2-index spinor, Z , is an 
eigenspinor ol* X , with zero eigenvalue, since 

= O (4.16) 

In addition, it can be shown that any antisymmetric 2-index spinor can be 

fl B 

written as a scalar multiple of the skew metric spinor, £ . Thus we 

need only consider the symmetric eigenspinors ''tj ^ of "% , that are 

then at most three in number and satisfy: 



Then since y and % are both symmetric in their spinor indices, the 
4x4 matrix system for the eigenspinors of % can easily be written as an 
equivalent 3x5 system. For if we define the new variables and X^ by 
the relations: 

^ \ , Y3 * %x 

and ~ , ^12. = Xrj » JiF ^mt j s X 

X/3 s X31 * Xiiii } X23 " X s J? 7 Cj ixi ^ X 33 = 

then (4.17) takes the form 

V j = •> y • < 4 - ,9) 

Thus the eigenvectors of the 3x3 matrix X are simply related to and in 1-1 
correspondence with the symmetric eigenspinors of the 4-index spinor % . 

Mote that X, like P, is a symmetric matrix. Equations (4.17) and (4.19) may 
be related by the index mapping: 
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A £ 


l I 


12, 2 f 

a 


2 . 2 . 

3 


and we may then write (4.15) in the form: 



l\ x • Z j 


( 4 . 20 ) 


( 4 . 21 ) 


where now all indices run from 1 to 5 and both P and X are traceless symm- 
etric complex matrices. (Note: X is symmetric only if all indices are up 
or all are down) . 


We now show that equation (4.21 ) actually represents a similarity trans- 
formation between P and X. Thus either P or X may be used in performing the 
Petrov classification. Equally, because of the relation between X and , 
the Petrov classification can be done in terms of the symmetric eigenspinors 
and eigenvalues of *% . 


The condition for (4.21 ) to be a similarity transformation is: 

A 


yA y* ^ q A 

4 i Z j *- ° f 


(4.22) 


This is most easily checked by again using the local Lorentz frame, in which 
the (S’ -matrices that enter the definition of £ can be chosen as mult- 
iples of the Pauli spin matrices. Note that, since T and B range only 
over the values 1,2 and 3> the index 0 is not used in the set ^UV|0 ^ . 

Thus we need only the forms of 6", (T* and 6"^ , for which we will use 

the representations: 


cl = k 


O I 


{% 



(4.23) 
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Using the definitions (4.15) we have: 



(4.24) 


Consider one term of this product. Raising and lowering the spinor indices 
using the metric spinor allows us to write: 


(fif sf-x 


Trace 


(4.25) 


In (4.25) we have used the relation $r ?Q s and we note that since 

the trace is taken in spinor space, the tensor index positions can be either 
up or down. To evaluate the trace, we use the fact that with the represent- 
ations of the 5 * - matrices used in (4.23), ~ . Thus we have: 


Trace (ff f £ C % £ (\ Ctt ) = Trace 

(45) 

The Pauli spin matrices satisfy the well-known relation: 

\ v + l /@ Cy 

Using (4.27), we find 


(4.26) 


(4.27) 


Trace (V) = Jfc ( - S/S /* + ?/ C) (4 ’ 2S) 

Performing the same calculation for each term of (4.24), and combining the 



results, we find: 


1H. 


= -S/S/) = S"s (4.29) 

r*' A" 

Note that since we require that j j = -1 when A = ( , V ) and 

B = ( V^^t. ), (4.29) is the appropriate form to use for £ £ . Equat- 

ion (4.29) confirms that the relation (4.15) is indeed a similarity trans- 


formation. 
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5. The eigen spinors of X- and the Debever-Penrose directions . 

7/e have shown in the preceding two chapters that the Petrov classific- 
ation can be described in terns of the properties of a complex traceless 
symmetric matrix P, and that this classification is invariant under a 
similarity transformation on P. Further, we have shown that P is related 
by a similarity transformation to the matrix X, and that the eigenvalues 
and eigenvectors of X are in 1-1 correspondence with the symmetric eigenspin- 
ors and associated eigenvalues of the wholly symmetric 4-index spinor X . 

It now remains only to determine explicitly the symmetric eigenspinors and 
eigenvalues of X and show how these relate to the Debever-Penrose directions, 
in order to complete the relationship of the Petrov method of classification 
to the geometric quantities of the Debever-Penrose principal null directions. 

First, since 7C ABCD is symmetric in all its indices, there exists a un- 
ique (except for scale factors) decomposition of 7C into a symmetrized 
product of 1 -index spinors. This result is an immediate consequence of the 
fundamental theorem of algebra, which tells us that a general quartic form 
has a unique factorization into a product of linear factors^ Thus we 
can write: 

'XrtBo * -£(* "ft Sj>) (5.1) 

Each single index spinor such as defines a vector in Riemannian space, 

from (4.l). This vector, say , is also readily shown to be null. Thus, 
the spinor 10 will in general define a set of up to four such null vectors. 
These are the Debever-Penrose principal null vectors, and to relate them to 
the Petrov classification we need to obtain the symmetric eigenspinors of % 
in terms of the single index spinors given in the decomposition of (5.1 ). 


20 . 


To find these eigenspinors, let us consider the product ■$* 

Writing out in terms of symmetries over just two spinor indices, this 
gives : 

^ filet = [-fciA-mj) + (c S j) ■+ 

+ \ h) 1n fe' f 'j>) + -t w m,) 4( C ?,) + SuW^+j) 

+ S,j -A fc -£( (5 ‘ 2) 


Noting that As = 4,s 2 -4iS, - -s„ 4 a 
and similarly for the other spinor pairs, we will write as 

and contract on C and I) in (5.2) to give us: 


(5.5) 

As 


'fcfiUcx A 

/ * 1V 51 "£(A f) { ^ + 

+ 't's) Sa-tw/ + 'ii(ASg) 4a 

+ ^' rt I) ^A-^. 4 a s + -vf A 4 

+ 'fy S$) ^*4 -£ A /yw (5.4) 


To simplify (5.4), we use the spinor identity: 

£ 7 s a>3 ) ^ 0 (5.5) 

? 

This indicates that the sum over all cyclic permutations of the 1 -index 
spinors k,m,r and s, symmetrized and anti-symmetrized as indicated, is 
identically zero. It is readily proved by symmetry arguments, but does 
not appear to be given in the standard reference works on spinor properties. 
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In particular, taking C = 1 and D = 2 in (5.5) we have: 


^ T ( ^B) ^A^) 55 Q ^.6) 

P 

To make use of this result, we regroup the terms of (5.4) to the form: 

7~/\6cj> ~ -$( A wp) ( «r A -& s A ^ + 

+ "tyl s 3) '*V A*4 + ^a 4 (-6^ -fj) S'aVv- -t S(A 1n 'B) “^Vf") 


and now apply (5.6) with the symbols ordered as k,r,s,m to the third term 
of (5.7 and ordered as k,s,r,m to the fourth term. This gives: 


X 


flScx 


n(c 1) 

'1C w 


( "^A "$> S / Jvv + -h -irv S'a -^ ) 

*(" ^ ^£A (^A^f 


Similarly, we find: 

%A 6 c.j> fV’ = + M S i)(^s A -n^ 4 'h.vu 

+ 3w5 M mj) (S ft t-T A S) 


Taking linear combinations of (5.8) and (5.9) then at once gives us the 
eigenspinors of lC ABCJ) : 






( 5 . 10 ) 


with associated eigenvalues 

"X" ~ <2 ('i'A S A "iru -+ Sa^) 

+ 6 (-^A^ s) 


(5.11) 



or 


we 


t 0 (c 3) 

Since we could just as well have started with Jfz ^ 
have as the full set of eigenspinors : 


^ in (5.2), 


1+ 

— S' ± s i) 


V 

= ~£(a^s) ^aS ± nm-iASg) 




(5.12) 


with their corresponding eigenvalues: 

^l + - + "/“a IK S A $) + (j 

'X l + - -ks^+ Tn A -rs A ^)? L{-k^i* h s) 

V = ^(^Tn^S + T 6(& a S ^ A 'tw) (5>13) 


The form of (5.12) suggests that there are six symmetric eigenspinors, but 

•j* I — M 

if we apply (5.6) we find that \ = ^3 » T/| = ^ . and 

. f — 

= / i/ . Thus there are at most three independent symmetric eigen- 

spinors, and we choose these to be y , = / y, + * = 'jJ’ » 611(1 

% - Vi , with the corresponding eigenvalues Tv, , 0^1 and . 


If now k,m,r and s are all distinct 1-index spinors, then in fact we have 
four distinct Debever-Penrose principal null directions, three independent 
symmetric eigenspinors of r %> » and correspondingly three linear elementary 

divisors of P, so space-time is Petrov Type I. However, it is possible for 
k,m,r and s to coincide in various ways and we need to examine the behavior 
of the symmetric eigenspinors (5.12) and their eigenvalues (5.13) in such 



. The possible situations that can arise are the following: 


k,n,r and s are all distinct. We have four independent Debever- 
Penrose directions; the symmetric eigenspinors of X and their 
corresponding eigenvalues are all distinct; P has thus three 
linear elementary divisors and three distinct eigenvalues, and 
space- time is Petrov Type I . 


If k = r ^ i / s, In this case, just two Debever-Penrose 
directions coincide. We have k^r = 0, so the symmetric eigen- 
spinors are: 

'TM'J) 'T'aS + ^(A Sj)| A T^ 

% - -ft = -*«■-**) m,s (5 * ,4) 

with eigenvalues = 4 “X, • 

Thus we have two distinct eigenvalues and two symmetric eigenspin- 
ors, P has one linear elementary divisor and two distinct eigenval- 
ues, and space-time is therefore Petrov Tyne II . 


If k = r = m ^ 
directions coincide. 


s, in this case three Debever- Penrose 


We find that 


\ • and \- 


x = x = 0. Thus we have one symmetric eigenspinor, P has 

no linear elementary divisors, and space-time is Petrov Type III . 


If k = m j! r = s, then the Debever-Penrose directions coincide 
in pairs. It is necessary to be rather careful in applying (5.12) 
here, since the form of the eigenspinors becomes degenerate. It 
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is thon in fact easier to re-derive the eigenspinors directly from 
(5.8), and one finds the three symmetric eigenspinors: 


1/j. ' $ 2 ) % = 



Thus % has three distinct eigenspinors but only two distinct 
eigenvalues, P therefore has linear elementary divisors and two 
distinct eigenvalues, so space-time is Petrov Type I-D. 


e) Ifk=m=r=s f all four Debever-Penrose directions 
coincide. In this case, we find the symmetric eigenspinors: 

^ jj) (5.16) 

where j is any 1 -index spinor independent of k^. The eigenvalues 
in this case are zero, so P will have one linear elementary divisor 
and zero eigenvalues, which means that space- time is Petrov Type II-N . 


For completeness, one final special case should be added. If % is identic- 
ally zero, then space is flat and we have no preferred directions. Anything 
i r then an eigenspinor, and the classification ceases to have significance. 


Although historically Petrov developed a formulation in which Type I-D is 
a special case of Type I, since both have linear elementary divisors, in 
terms of the Debever-Penrose directions Type I-D appears much more like a 
special case of Type II, in which a second pair of Debever-Penrose directions 
are allowed to move into coincidence. 
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This completes the general development in which we relate alternative ways 
of classifying space-times. For the remainder of this work, we will conc- 
entrate attention on algebraically special space-times, and utilize the 
structure that this implies to look for solutions of the field equations. 


One final point should be made here. For ease of development, we have 
worked always with the vacuum case. When matter is present, the Weyl 
tensor is used in place of the Hiemann tensor in the Petrov classification, 
and a corresponding gravitational spinor, 'ftscJD , is used for the spinor 
formulation^^. Then the analysis goes through in a way that exactly 
parallels the vacuum treatment. However, in terms of the utility of the 
results, the use of algebraic degeneracy has proved fruitful in looking 
for solutions mainly in the vacuum case. Although solutions of the 
Kinstein-Kaxwell equations have been found for algebraically special 
space-times^’ no solutions of great physical interest seem to have 
been discovered for such situations. We will discuss the non-vacuum case 
further in later chapters, and point out some of the added factors that 
complicate the solution of the field equations in such cases. 
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6. A summary discussion of degenerate space-times . 

In any space-time which is not Petrov Type I, two or more Debever-Penrose 
principal null directions must coincide. Since this is a statement about the 
essential geometry of the space-time, it is natural to ask how this geometric 
statement relates to other geometric properties - for example, must an algebra- 
ically special space-time possess some physical symmetries (isometries)? Counter- 
examples show that physical isometries neither imply nor are implied by algebraic 
specialization of a space-time: algebraically special spaces exist that have no 

physical symmetries, and there are algebraically general spaces that possess 

( 35 ) 

one or more isometries . 

The best insight into the meaning of algebraic specialization is gained by exam- 
ining analogies with the electromagnetic field. However, the nonlinear nature 
of the gravitational field equations makes many of the methods often used in 
handling the Maxwell equations ( such as Fourier decomposition) inappropriate^ . 
To summarize a substantial body of work in a few sentences, it turns out that the 
most enlightening analogy is in terms of the eigenvectors (eigenbivectors, in the 
case of the Riemann tensor) of the field tensor. In the case of the Maxwell field, 
only two situations occur: the eigenvectors of the field tensor are independent, 
or they coincide. The former case occurs in any region with radiative sources 
present, the latter describes a pure radiation field asymptotically far from 
bounded sources (34*36,37)^ 

The gravitational field tensor has a more elaborate possible structure for the 
eigenbivectors. Five cases can occur, depending how the four Debever-Penrose 
principal null vectors coincide. Thinking in terms of gravitational radiation 
from bounded sources, five different possible cases (which are related directly 



to their Petrov types) multiply different powers of a mean inverse distance 
from the bounded sources. This "peeling-off" property is of interest 
to us here only because it indicates that gravitational radiation in regions 
containing sources cannot be Petrov Type II or Type III - i.e. snace-time near 
radiating sources must be algebraically general . 

This result suggests that solutions obtained for algebraically special space- 
times will either be non-radiating ( like the Kerr and Schwarzschild solutions) 
or applicable only to pure gravitational radiation far from all sources. It is 
interesting to note that the Kerr and Schwarzschild solutions are of Type I-D, 
which does not occur in the "peeling-off" theorems, and seems to be associated 
with stationary solutions rather than with gravitational radiation. Very recently, 
other stationary asymptotically flat solutions that have non-zero angular momentum 
and are Petrov Type I have been developed, but their possible physical signif- 
icance is not yet clear (39»40)^ 

Despite the fact that algebraically special space-times cannot describe the most 
general physical situation, the problem is still sufficiently general that as yet 
no one has succeeded in constructing the most general form for the metric of an 
algebraically special space-time. It is quite easy to define a necessary and 
sufficient Condition in terms of the Riemann tensor: the matrix P of Section 3» 
Equation 10 , must have at least two equal eigenvalues. Thus if the character- 

istic polynomial of the traceless matrix P is: 

(*x) ~ + c (6.1) 


then space-time is algebraically special if and only if t(A) = 0 has two equal 


root3, which requires that: 


27 c 2 = 4 b 3 (6.2) 

This elementary and apparently attractive relation between b and c becomes an 
extremely complicated nonlinear partial differential equation in the metric 
tensor when we substitute the forms of b and c. As a result, the direct approach 
of (6.2) is quite useless. Instead, most efforts have been devoted to the analysis 
of particular classes of metrics that can be shown to correspond to particular 
algebraically special space-times, without seeking or claiming the most general 
possible form. 

In the subsequent chapters, we develop a very elementary approach to just such 
a class of metrics, first studied using a tetrad formalism by Kerr and Schild^ \ 
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7. A class of metrics with algebraically special space-times . 


We now focus our attention on a particular class of space- times, with metrics 


of the form: 




(7.1) 


where in (7.l) is the Lorentz metric, m is any real arbitrary constant, 

and styi is a 4-vector that is null with respect to the Lorentz metric, so 

* 0 . 

It can be shown directly^ ^ that the metric (7.1 ) corresponds to an alg- 
ebraically special space- time, but we prefer to establish this in the course 
of our general development. 


A few comments are in order on the form of the metrics we are considering. 
First, if we make the choice 

^ = (V) * (l, > Z /r) (7.2) 

2- « • o 

where ^ = 04 .+• 'vj -*-2 , then we obtain a line element that is the 

(42) 

Eddinington form of the Schwarzschild metric, obtained from the original 
form of the Schwarzschild metric by a change in the time coordinate. In this 
case, the arbitrary constant m appears as the mass of the body. Second, if we 
are interested in asymptotically flat space- times , must tend to zero at 
spatial infinity if are Cartesian coordinates for Minkowski space- time. 
Third, for small m and finite , (7.1 ) has the form of a perturbation on 
Minkowski space, thus we expect that the usual linearized theory will emerge 
for small m. 


In the treatment given hereafter, we will treat the spatial coordinates in a 
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symmetrical way. We will also make no assumption of stationary or axisymmetric 
solutions, nor do we introduce the assumption of algebraic degeneracy explic- 
itly in our analysis. This differs from the treatment of Debney, Kerr and 
(43) 

Schild , who in a tetrad formalism use the algebraic degeneracy from the 
outset, and from that of Misra^^, who assumes axial symmetry. 


Using (7.l), a number of useful relations are readily established. We find 


g/* v = -Tj / 4 ' 1 + a, 


(7.3) 


so that indices on may be raised or lowered with either the full metric 

tensor or with the Lorentz metric tensor. Also, since is null 




and from the product rule for covariant derivatives we also have 




'/Lit- 


j 


= 0 


(7.4) 


(7.5) 


The Christoffel symbols are easily calculated and we find 


(7.6) 


The field equations are greatly simplified by the choice of metric, since 


from (7.l) we find det(g) = -1, and thus we have 


lL] » (WA 


- 0 


(7.7) 


Using (7.7), we find that the field equations reduce to only two terms 


7v<v 




- 0 


(7.8) 
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The metric tensor is linear in the arbitrary constant m, thus is a 

fourth order polynomial in m. If we require that (7.1 ) shall lead to a sol- 
ution of the free space field equations for any value of m, then each power 
of m in (7.0) must vanish separately, which gives rise to the four sets of 
equations: 
order m: 

= 0 (7 - 9> 

order : 

3 

order m : 


4 

order m ; 

= O (7.12) 

Using (7.3) - (7.5), we can at once verify that the order equations are 

3 

satisfied identically. The order m equations, after expansion of the 
Chri3toffel symbols, can be written in the form 




o< 


(7.13) 


where we have defined /\F S JL JL latter is null, and it is 

also readily seen to be orthogonal to the null vector -d. . This can happen 

only if / VT and are proportional to each other at every point, thus 

w'e can write 


ar~ = JirjTu = - U' 


(7.14) 
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where A is a scalar field. 


Note that (7.14) tells us that the vector field is tangent to a family 
of geodesics' , and is parallel-propagated along the geodesics. If, 
as will happen later, we choose to perform a scalar change of variables 
such that (7.14) becomes 


oi 


\y< * -re.' i 


- « 0 


T 


(7.15) 


then we say that the vector field is affinely parametrized. There 

r 

are many possible choices of the scalar H that will affinely parametrize 
a given family of geodesic tangents. 


Now considering the order m equations, we define a new scalar L by: 


L = -£ 


jo< 


(7.16) 


and then we can rewrite (7.9) in the convenient form: 




i'Lu-C-j) 


Ik 


L ^z. 


(7.17) 


f— I ^ ^ r-** 

where the D’ Alember tian operator is defined as \_J = /j£. l ""V » s0 

the upper index derivative in (7.17) is to be raised using . 

Defining G = (L + A ) (7.18) 

we have as the final form for the order m equations: 


- (Vev) V = 'A (<?-ty), y) 

These will be used extensively in subsequent chapters. 


(7.19) 


For the order m equations, expanding the Christoff el symbols and using 
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equation (7.14), equation (7.10) can be rewritten as 


2 + 4, " t — A 2 - 0 (7.20) 

Using the definitions of L and A, we readily derive 

^\f = [( L ^ (7.21) 

and by using (7.19) to form we also find 

~ (L Z 'A 1 ) + ((P^/jj^ (7.22) 


Substituting (7.21 ) and (7.22) in (7.20), we find that this gives an 

2 

identity, thus the order m equations contain no new information and 
are implied by the order m equations. This is expected on general grounds, 
since there are ten independent order m equations, and this should suffice 
to determine the metric completely. In the same way, we can also show that 
the relation (7.13) leading to the definition of also follows from the 

3 

order m equations, although the derivation from the order m equations is 
shorter and simpler. 

We now work exclusively on the order m equations. From (7.19) 


p ( no sum on ^ ) (7.23a) 

- kv «,r. M = a. ( c«v)iv (no sum on y ) (7.23b) 

using (7.23), we can now eliminate the second order terms from (7.19), 
which leads to the six independent first order equations: 
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y^l^^y ^u/c^-^v^^y 2, — -£y|o< ^y |a< ^f 

^ {A'^aIv) ^ -Sy 2 — -^V|v-^/vJ 


(no sun on and V ) 


(7.24) 


(7.26) 


These first order equations, together with (7.23a), are completely equiv- 
alent to (7.19). If we now take 

(7.25) 

where H is any scalar field, we find the form: 

-jj- (y*|v) 

(no sum on and v ) 

This has exactly the same form as (7.24), but bas replaced A /" 

and g/h has replaced G. We also readily confirm that is null and that 
(7.4) and (7.5) are satisfied with replacing -^u . If now in (7.26) 

we make the particular choice H = st 0 , so that -& 0 = 1 , then we have: 

■ 2 -«- 4 x-Hi'Ui-kj - 4;u4i u -t/ -4 jU 4^4- 


(no sum on i and j) 


(7.27a) 




(7.27b) 


(no sum on j) 
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Roman indices in (7.27) and subsequently run from 1 to 3. 


Using (7.27b) to simplify (7.27a), we find a form completely equivalent 
to equations (7.27): 


- l * - -JE 

4 




(7.28) 


Finally, using (7.25) with H = -d 0 in (7.14), we have the relations: 

= — A (7.29a) 

- -(U'4* (7.29b( 

Eliminating A between (7.29a) and (7.29b) at once gives: 




\y. 


= -kt'-ffy = 0 


(7.30) 


Thus is an affinely parametrized null vector. We note that the choice 

H = st , where JL is any component of M , also leads to an affine 
parametrization. However, the choice of the time component is the most 
natural, since it permits a subsequent analysis that is symmetric in its 
treatment of the spatial coordinates. 
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8. The reduction of the first order field equations . 

Equations (7.23^) and (7.28) comprise four second order and six first 
order differential equations, which together are equivalent to the 
entire set (7.19) of second order field equations. To solve the field 
equations, we first address the reduction of (7.28) to a simpler form. 
Noting that we can lower indices on both and JkT using , we can 
write 4 = — . We further define p = ^/%t 0 > so that 

equations (7.28) become:- 

-$ulo-£j|o — 'fcu-o -4jU = 

-t>Oii ♦ *iii -[4^4*] (8 -’ ) 

where we are now summing over all repeated indices. 


Noting that (7.30) can be written as: 

•"^lo - £iu H 14 , 

we write (8.1 ) as: 

- |> [Ui-ftj)i4 4 t - + 4jic)'] 

If we now define a 3x3 matrix M and a 3- vector ~fk by: 

Mtj - -An j 

fit); = 4i 


( 8 . 2 ) 


(8.3) 


(8.4a) 

(8.4b) 


then (8.3) can be written as the matrix equation: 
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'/j,[n( i - -iU T ) n T ] = 

m+m t - [n£.)£ T - £(n£) r (s - 5 

We now perform a rotation of the coordinate axes, so that 'fi goes to the 
form Jfc. - I 1 \ . 


In these coordinates, using the relation 

n , t 4' = n T £ -4*4, - o 

I 

we find that M must have the form: 


h 


/T 


0 ' 

L - -i 

0 I 
I 

I 


'V 


V 


m* 


where m’ is a 2x2 submatrix. Hence: 


( 8 . 6 ) 


(8.7) 



Using (8.7) and (8.8) to form (M' + K ,T ) and M'lc' (M , £ , ) T , 

and substituting in (8.5), we find that x and. y cancel from the equation, 
and we are left with the remarkably simple result: 


On‘W T = J? f/yn' + W T ) ( 8 . 9 ) 

Equation (8.9) implies that m' can be written in terms of a real 2x2 
unitary matrix U, as: 


U 


X - *J/,, 


(8.10) 
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For our main development, we will assume that U is proper, with positive 
determinant, and so can be written in terms of a single real parameter 
0 in the form 

1 CoS 6 -SlrvO 

y CM 0 

The cases of improper U, vanishing p and vanishing 0 are discussed in 
Chapter 13. 




Using (8.11 ), M' can then be written as: 

[) I 0 0 0 \ 

U \- CmB "Svw0 

V S^v9 j 


= 




( 8 . 12 ) 


where we define pu = x and pv = y. 

To get back to the original coordinate system, we must form 

M = RVR 


( 8 . 13 ) 


where R is a real 3*3 orthogonal matrix. Using the orthonorraality of the 
rows and columns of R, and noting in particular that 


we find the form for M: 


( 8 . 14 ) 


( 8 . 15 ) 



Now defining 


T; - \> (u X u + v^ 3 c) 


(8.16) 


we note that r.ince and R^ are row3 of H and therefore orthogonal 
to R is orthogonal to K^. Further, using = 0, we 

also have 


~Ti -Al = 0 

— k 

so T is orthogonal to k. 


(8.17) 


Using (8.17) in (8. 15) and multiplying by kjc.., we find that k^R^ = 1. 
Since both 1c and ”R^ are unit vectors, we must have *k = + T? , and we 
can choose the plus sign since the overall sign of R is arbitrary. 

Then using (8.4a) we can write (8.15) as: 


M = «*(V-M,-)+ 

where in (8.18) we have defined and jl by: 

ok' — j> ( J “* o>& &) 
p = f> S^9 

Multiplying by k. and using (8.2) we have at once 

~n = 4;,. 

so (8.18) can be written as 


(8.18) 


(8.19a) 

(8.19b) 


( 8 . 20 ) 


+ ~&l\o -&j 
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( 8 . 22 ) 


or, using (8.2), in the alternative form: 

-&U ( - -^j) = °<r l 

The forms (8.21 ) and (0.22) will be used extensively. They replace the 

six nonlinear first order equations of (7.28) by an explicit expression 

for the derivatives k. . . Before studying these equations further, which 

1 1 0 

is our task in the next chapter, we first note an important feature of 

(8.22) . The matrix ( £ 1 ; ~~ ) is singular, thus although 

vj ^ J 

(8.22) is a general set of relations between the k. . . and the vector 

1 1 J 

k., we cannot use (8.22) to solve for k. , . in terms of k. . The only 
i i|j i 

exception to this statement occurs if k. . . k . = 0 , when (8.22) becomes 

soluble for k. ( .. However, from (8.2) we see that this occurs only when 
i[ J 

^i|0 = which is the stationary case. 


For the general non-stationary case, the additional information needed to 

solve (8.22) for k.. . will come from the equations (7.23a) which we have 
i| J 

not yet considered. 
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9. Properties of the first order equations and the optical scalars . 

The equations (8.21 ) lead to a remarkable number of interesting relations 
involving and S , but we will confine our attention to the main 


results needed for subsequent analysis. 

Using k^ k_^Q = 0, (8.21 ) at once gives 

= (9.i) 

Also we find, directly from (8.21 ), the following:- 

£«,je -Jl ^|j — <2 ft (9.2) 

-&ilj- 6 jK ~ ( 9 . 4 ) 


In the study of geometrical optics in a gravitational field, the quantities 
known as the optical scalars play a fundamental role^'*'^*'^. Two of 
these scalars, the expansion 6 and the twist QJ » are defined in terms 
of an affinely parametrized ray vector field by^:- 

6 = % 4 r \y> (9.6) 

w = [/ 3L Or^iw - £vij^4/ 4W ] 4 (9 - 7) 

Using (9.1 ), we readily show that: 


(9.8) 



and also that: 


Now using (9.3) and (9.4), we have: 

-4vf^04^ ,v = ‘Sj? 2 

Comparison of (9*6) and (9*8), and of (9*7) and (9*9:) indicates that except 
for a choice of sign and are exactly the optical scalars 9 and 

D . It is interesting to see that the optical scalars enter the solution 
in a natural and fundamental way, without being introduced at the outset. 

The third optical scalar is the shear, 6* . For metrics of the form (7.1 ), 

the Goldberg-Sachs theorem tells us at once that the shear must be zero, if 
the space-time represented by (7.1 ) is algebraically special. Conversely, 
a direct calculation of 6T for the metric of the form (7.1 ) yields 6 s = 0, 

which confirms directly that these metrics correspond to algebraically special 
space-times. 

(The shear is defined in terms of an affinely parametrized ray congruence 

. . . (34) 

k / to be : 

lfrl = ( 9-"> >• 

Using (8.2), (9.l) and (9.4) to calculate leads to an important rel- 

ation for oc 

^iv « (W**- p X ) 



(9.9) 


( 9 . 10 ) 


( 9 . 12 ) 



Similarly, using (8.2) and (9.2) and calculating and 

leads to an analogous relation for |3 

fly V ~ <2 oSp 

Introducing the variable 

X =• + if 

(9.12) and (9.13) combine to the single equation: 

V 




(9.13) 


(9.14) 


(9.15) 


The complex variable Y plays a fundamental role in the following disc- 
ussion, and (9.15) will be used frequently. 


One of major objectives in the subsequent development will be to derive 
differential relations involving k^ and Y . However, there is an alter- 
native method of seeking a solution, which directly explores the integrab- 
ility conditions of (8.21 ). We will not discuss this method in detail, 
since it forces us to abandon a symmetric treatment for the spatial 
coordinate variables, but in Appendix 1 we give a brief discussion of the 
approach for the stationary case when k^j^ = 0 . 

To develop differential relations, the 3-gradients of and |3 also 

prove useful. To obtain, these, we use (8.21 ) to form k, . .. . , and we 

^ I J I J 

also form )| s using (9.5), to give us a second equation 

involving . Eliminating j j | j between these equations, using 

(9.l) and applying (8.21 ) to evaluate k^j^ ^j|i we the form: 
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+ ZfU> t 4 e 4it, + % p,j^ t (9 - ,6) 

Forming £ t from (9.16) and using (9.15) then leads to: 

fit ~ ~c2otp-£t -*• 2f-fUio ~ £ 4 ^ ^{j ~&f> ^ 9 * 17 ^ 

Equations (9.16) and (9.17) cannot at once be combined to a simple equation 
involving only , and we will later show that we must introduce a slight- 
ly different Y to achieve such an equation for the gradient. For the 
moment, we simply note that it is only the terms of (9.16) and (9.17) 
involving k^j ^ that fail to combine to a simple form in ^ , and thus 
an analysis of the stationary case should be possible without redefinition 
of this complex scalar. This is precisely the path that we followed in 
treating the stationary case in Reference 27. 



10. Reduction of the second order field equations . 


We now turn our attention to the manipulation of the second order field 
equations (7.23a). We look to these equations to provide the information 
we need to resolve the indeterminacy of (8,21 ) and (8.22), thus we expect 
that some relation between the time and space derivatives of k^ should 
be provided from the second order equations. Using = J@ 0 , 

equations (7.23a) become: 

!* = A(frOi. 


U/), 


(10.1a) 

(10.1b) 


(no sum on i) 

If we add the three equations of (lO.lb) and subtract (lO.la) we have after 
using (9.1 ) : 


( (r-'do)^ = [G-@d] (10.2) 

Consistent with our comment following equation (8.11 ), we may now assume 
that ^ 0 . Comparing (10.2) and (9* 13)# it is natural to take 

Gz 0 - Xp (10.3) 

which at once leads to the relation: 

r= 0 (10-4) 

\ p /(yH. 

If we now take (7.27b:) and sum over j, we have: 



- 


(10.5) 



Using (9.5) now at once gives 


Q = -k[^\e z ) 

* 

Now using (1O.4) and (10,6) together with (9.12) and (9.13) gives: 



This suggests that we should write, defining W : 


W ex' - £0 


( 10 . 8 ) 


and then 


VI y> 4 ^ - o 

where W is a real scalar variable. 


(10.9) 


Making use of ( 10.6) and (l0.8), the field equations (7.19) take on the 


form: 


+ [ (10.10) 

Taking 0 then leads to the four equations 

(10.11a) 


Using (7.28), it is then straightforward to show that ( 1 0. 1 ) and (lO.1l) 
can each be obtained from the other and we will therefore work with (lO.1l) 
since they have the advantage of being linear in k^ . Subtracting (l 0.1 la) 
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times from ( 10. 11b) then gives us the equations: 

+ a Utp -fh'iV + s.w^4 1 1 ! 1 

= -[w^+jn],. 

- W (oS l + o 


( 10 . 12 ) 


We now evaluate each term of the left hand side. Using (8.21 ), (9.12), 
(9.5) and (9.17) we find: 


^1(0-^^/° - \ ^*)f£ + 


lo 


-H 


■^ujo ( J3 1 ) - (u\ f l ) 


(10.15) 


Similarly, using ( 1 0.9) leads at once to: 


IV| 


— ^ Wji, ■+ 4i\rJ l0 - p W|j ^ 

(10.14) 


and finally by differentiating (8.21 ), and then using (9.12), (9.5) and 
(9.16), we find: 




(10.15) 


Using (10.15), (10.14) and (10.15) in ( 1 0. 1 2 ) at once gives: 

(p - ^ ) ( \7jo -fee - lV(i,) + 3-^f> £^1 V!\^4ji 


(10.16) 
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or, putting W = P - , 


= -'P(y‘+p‘J-€i|. 


These are three first order linear equations in P. We will use solutions 
of these equations to provide the additional information needed in (8.22) 
to determine k. . . in terms of k. . 

In particular, we note that if P is a function of alone, and has no 
explicit dependence on the coordinates, then (l0.9) is automatically 
satisfied since 

W 21 (4;u4. / i = 0 (10 ' 18) 

cL'p * 

With the assumption that P = we can simplify ( 1 0. 1 7) using 

(8.21 ) to the form: 




( 10 . 19 ) 


Consider now the possible forms for P. It is to be a real scalar function, 
depending only on . We cannot construct a scalar by contracting 

on itself, since is null, and if we contract with any other 

4-vector, the latter must be independent of position. This strongly sugg- 
ests that P must be a function in which hL is contracted with a set of 


constant 4-vectors, thus: 




( 10 . 20 ) 
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where 


a 


is a vector with constant components. Substitution of any such 


form as (10.20) into (10.19) will lead to a relation between k^ r and k^ ^ 
We now 3eek the simplest such relation, by supposing that P can be written 
as a Laurent expansion in terms of a single scalar , thus : 


T ■" ^ ( /fxA 


(10.21) 


ov s -oo 


Each term in the expansion, used in (10.19), leads to the form: 


-' nU^Hs = [ a .\- k * Ax ) A llo O0.22) 

Thus any value of n gives a relation between ^j_| r ^ijO ’ we see 

from (10.22) that if and only if n = +1 , we obtain a relation that does 

not involve k itself. Taking n = +1 gives the very simple result: 
r 


or 


~ & -&LU 

/*£ = 0 


(10.23) 

(10.24) 


We then have for P and W : 


? = 


T 


W = [*r-& r ) 3 


(10.25) 

(10.26) 


r 

Since P is real, it follows that also has real components. 


We will shortly return to consider (8.21 ) for the particular form of P 
given by (10.26), but before doing so we make one important observation. 
In order for a space to admit a Killing vector /erf* with constant comp- 
onents, we must have 
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(10.27) 


For the metric we are using, from ( 1 0.8) we see that (10.27) will be 


satisfied if: 


/«7 — 0 ( 10 .: 

Using (9.1 )» (9.2) and (10.2?) we readily derive the following: 

Wu saf ■= O ( ,0 - : 

= o (10 - : 

Thus (lO. 28a) is at once seen to be satisfied, and using this and (10.23) 
in (lO. 28b) we see that (l0.28b) is satisfied. 


This implies that /&/*’ is a Killing vector of both and of the back- 
ground metric * a fact that we will use later to simplify the final 

equations for • 
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1 1 . Defining differential relations for the function Y . 


We are now in a position to establish the basic differential equations 
that the complex scalar Y must satisfy. Using ( 10.23) in (8.21 ) 
gives us: 

/Ot -$t i, jj ■+■ “ ^ ^ j 

~ (ii.i) 

0 , 

Assuming for the moment that (CK, $ 0 , an assumption that we will 
examine later, equations ( 1 1 . 1 ) become: 


{'CL + yCL^i.^) =r 

/CL (£ij Hi j ■+ /CL p CujJt -dj 


( 11 . 2 ) 


The solution of (11.2) for r is simple, since the inverse of the matrix 
(/ ) is O ust ( — 'ctyj /f '01? » readily verified 
by direct multiplication. 


Thus (11.2) is equivalent to: 

= [W (£u— Ac- 8 /) + 

■ X [? + j - -tf-ftf/p] 0, ‘ 5> 

As we remarked following (9.17), the scalars and |5 have the 

disadvantage that their gradients do not conveniently combine to a form 
in the complex scalar . We therefore find it convenient to introduce 

at this point rtew scalars ^ and R , which we will show have suitable 
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gradients thnt combine to a 
new variables by: 

CX? - 

? - ?f 


form in )f = ex' -t i . 


We define the 


(11.4) 


■h 


We will also define a lower index form of /Cl by the relation 


/C^m, = /Ct (11.5) 

We should note here that indices on /&/* cannot be raised and lowered 
in general using Ov ^ . However, we will later find that all the results 
we derive have exactly the form they would have were we to use 
instead of to raise and lower the indices on . As a pract- 

ical point, all raising and lowering operations with the metric of the 
form (7.l) seem to be accomplished just as well with as with 

which again emphasizes the very special form of this metric. 


Using (11.3), (11.4) and (11.5) then gives us: 

”^•'1 4 ~ Sij + P SyY-^e.) 

+ (£■- -g { 0°)^ + (,1.6) 

In reaching equations ( 1 1 . 6 ) , we assumed that /Cl ^ 0 . However, if 

sr 0 , direct calculation from (10.23) and (11.2) gives us exactly 
the same form as ( 1 1 .6) , as is shown in Appendix 2. Thus we can proceed 
to use (il.6) in all cases. 
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Equation (ll.fi) serves as a replacement for equation (B.2l), and will allow 
us to obtain Governing equations for V . We will make heavy use of it in 
what follows. 


Using (9.l), we now consider: 

S(T^) \jl - -^icnu 


(11.7) 


After a rather lengthy calculation, given in Appendix 3, we derive the 
result: 

- (^|0 + ^ ( f 

+ (^x|j -f- (4 x yk) 

r (11.8) 

where we have defined "a = ^ • 

Forming Vol x -{L * we at once find: 

Vp = (p|, - a. 2 3(3)4. + ' 2L*J 2. 

( »< ^ )( /GL X -$.) -h X (11.9) 

Using (9.12) and (9.13)» the last two equations can be written: 

V? ® (?,. 

i' ( Vy x ) — i y (/ctx'&j (11.10 


The above very useful result is the anticipated replacement for (9. 16) 
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and (9.17), and takes on a somewhat simpler form in terms of the new 
variable 



thus: 

Vw = ( (j O io -+ / a°)^i — ~X 

+ L(Vv x -tC) + (11.12) 

Now using the relations 

Vu> . = a 3 , 0 ^t° (11.13) 

Vu) • ^ = c 5| 0 + T (11.14) 


which are derived easily from (9.15) and ( 10 . 29 ), we find: 

(□u> f = I0| v CO ^ — ~ /OT (1115) 

The second defining relation satisfied by and U) is found by 

taking the divergence of (ll.io). Again, after a rather lengthy calculation 
given in Appendix 4 we find the simple result 

□ 2 y = y | V |V = 0 ("-'6) 

or in terms of U) 

U) D 2 W = — A ( /of 1 & uv) (11.17) 

Equations (11.15) and ( 1 1 .17) are the two fundamental equations satisfied by 
our complex scalar generating function. Before considering the equations 



in fact determine the metric 
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further, we will confirm that LO and ✓a/ 1 
completely. From the form of the metric, and from equations (10.8) and 
(10.26), it is clear that /CL and determine the metric completely, 

thus it is sufficient to show that u) and iGf completely determine ■hr. 


Writing (ll.ll) in the form: 


[Vu> + ) =• ( U 

and now writing 

ho -+ <L° + i- ( Vu> + & 

(11.18) 

— * 

Q 

^7 (a) H 47^ 

(ll. 19a) 

S 

- U )£> + sCL 

(11. 19b) 

we readily verify that 

Q * <5* = SQ*) -f 34 

— V — 

where B is a function of k and U) . To find B, we use 

(11.20) 


0 * (Q + Q*).(Qx cC) = 

sq r ).(«4Q*) + 3 (q + q i '').-4 ("-2') 


Prom ( 1 1 . 1 8 ) we f ind that 

(Q Hr GT*). ~ S* (11.22) 


thus we have at once that: 


3 


( S*Q, Q -+ SQ*Q*) + $.$*(StS*)] 

• 4 - (s +S*) 


(11.23) 
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and thus finally 

4 = l ( S $*+ s '$)]/ 3 ( 11 . 24 ) 

Since B, S and Q depend only on CO and of f we have confirmed that ( 1 1 • 24 ) 
determines At uniquely in terms of U) and /Ct/'. 
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12. Use of the Killing vectors and discussion of solutions . 

The equations for the complex generating function 0 can be simplified 
by use of the Killing vector At . There are only three cases to consider: 

a) If PC is timelike , then a Lorentz transformation must exist that takes 
it to the form (l,0,0,0). In this case, equations (10.26), ( 1 1 .15) 
and (11.17) reduce to the forms: 

W - 1 (12.1) 

(Vw ) 2 » 1 (12.2) 


w - X 


(12.3) 


Since there is now no dependence on the time, we would expect to 

obtain this case directly from (8.21 ) with set equal to zero. 

As we remarked earlier, in this case the first order equations are then 

explicitly soluble without using the second order field equations of 

/ \ 

(7.23a). This approach has been carried through in detail'' ' and 
does in fact yield the exact equations ( 1 2.2) and ( 1 2 . 3 ) . Note that, 
since P = 1 , the redefinition of the variable Y in this case 

is the identity transformation, and hence \ is the same as V • 

The most important known solutions of the stationary vacuum case are 
now obtained very easily. First, let us note that the general solution 

/ j ry \ 

Of (12.2) can be written using the theory of envelopes' in the form: 


LJ 


= Ai X; -+ b (Ai) 


(12.4) 
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where we define: 

A, 



(12.5) 


where b(A^) is any function of the , and where A^. A^ = 1 . 

The simplest possible choice we can make for b(A^) is to set it equal to 
zero. Then we have 


where 




r 


x ± . x . This gives us 


u) 



W oC = r _1 


( 12 . 6 ) 


r , so cxC = r and thus 


(12.7) 


We now calculate le immediately from (11.18), using = 1 . (Normally we 
would use (11.23) and (11.24). However, when u3 is real, we can see at 
once from ( 1 1 .18) that we must have 




-+■ /<*. 


( W (o + 


( 12 . 8 ) 


and thus in this case 

•I = vs = + U 

From (12.7) and (12.9) we then have 

-*r = fcf* ( 1 > % , % > % ) 


(12.9) ) 


( 12 . 10 ) 


and comparing this with (7.2) we see at once that we have produced the 
Schwarzschild solution with the Eddington form for the metric tensor. 
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By direct substitution, we also readily verify that the choice we have made 
for U in this case also satisfies the second equation ( 1 2 . 3 ) - 

A second simple choice for b(A^) is to make it depend on just one component 
of A, thus: 

( 12 . 11 ) 


( 12 . 12 ) 


(12.13) 

and again we readily verify that this choice satisfies ( 1 2.3) also. 

If K is a real quantity, then the solution corresponds to the Schwarzschild case, 
in which we have merely displaced the origin of coordinates along the z-axis. 
However, if K is imaginary, then use of (11.23) and (11.24) gives us a new 
result. Writing K = ia , we find: 

= (f* + -»tf)/u* + f) ) 

-fit = -■*■*■)/ { A. v +f> i ) > 

-4, * J 


(A; ) = 1C f\ 3 

where K is any constant. Using ( 1 2.5 ) then gives us 

A, = */ H t At - Vu 

A 3 = (z +t :)/w 

N = ( x *+ (z+fc) s ) !i 

Thus using ( 1 2.4) we find 

U> s iO * ( + ( Z--+ X) l ) X 


(12.14) 



60 


where yO is the real part of UJ and is a solution of the equation 


(18). 


(** ~~ f ~ — /Ct? Z* = 0 


(12.15) 


Similarly, we find 
2 




r7( 


i al 

fi + 42 j 


(12.16) 


and using ( 1 2. 1 4 ) and ( 1 2. 1 6) in (7.1 ), we find the form for the line 
element to be : 

2 i . r t i ,i_ , n o 3 


r — doc 1 — — — o f 

ft L 

+* + z^z/^ 


(12.17) 


This is exactly the form for the line element that was given in Kerr’s 
( 1 8 ) 

original paper' ' . A great deal of work has subsequently been done on 

both the geometry and the implied physics of this solution, and it would 

be inappropriate to reproduce much of that here. We merely remark that the 

form (12.17) is not the most revealing form for the line element, and for 

future reference we will quote an alternate form derived by Boyer and 
(48) 

Lindquist , in which the axial symmetry of. the solution is quite 
explicit: 


(Us 2 = jdt * - 4- A,Si^0 

( -f"*+ ^(^d) 

- ^+V(^‘-2^ + 2)) <l2 -' 8) 
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This form in obtained from the Kerr form by transformations of the time 
and the azimuthal angle variables^^’^*'’^ and we will find it very 
useful when we consider interior solutions in Chapter 14 . 

As is readily seen from ( 1 2 • 5 ) * most choices that can be made for b(A^) 
lead to nonlinear equations from which A^ must be determined, and to date 
no other solutions of physical interest have been found for the case with 
timelike Killing vector. 


b) If /at is spacelike , then a Lorentz transformation exists that takes 
it to the form (0,0,0,l). In this case, equations ( 10.26), ( 1 1 .15) and 
(l1.17) reduce to the forms: 

W - 


_ / 3 w 


r -a.T 


L 0 


VS - 


(SI -(If - 

i? — 1 s 

** 3^ J 


(12.19) 

(12.20) 

(12.21) 


This case has no dependence on z and thus represents a 2-dimensional 

problem spatially. Solutions of this form are not candidates for the 

representation of gravitating bodies, though they are appropriate to 

certain cases of gravitational radiation, and cylindrical radiation has 

( 51 ) 

been extensively studied ' . 


c) If &t is lightlike , then a Lorentz transformation exists that takes 
/Of to the form (l,0,0,l). We then have, from (10.26), ( 1 1 .15) and 

(11.17), the equations: 
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where 


W = d/U + ^f 





6J ;= 



= O 


W ( *, fc -*) 


( 12 . 22 ) 

(12.23) 

(12.24) 

(12.25) 


From (12,23) and ( 12.24), 10 must be a function of (x + iy) or (x - iy) . 
Using this in ( 1 1 . 23 ) and ( 1 1 .24) permits an explicit calculation of k^ , 
and leads to : 


k 1 = k 2 = 0 ; k^ = -1 


( 12 . 26 ) 


However, this means that W as defined by ( 12.22) becomes indeterminate. 
Thus we are led to conclude that there are no solutions of this form with 
a lightlike Killing vector. 


This is a surprising result, since previous papers by Kerr and Schild^^ 
and by Dabney, Kerr and Schild' ' give classes of solutions for metrics of 
the form (7.1 ), and these solutions include cases with lightlike Killing 
vectors. At first sight, the results given in these references are quite 
unrelated to the governing equations given here. However, one can show 
(see Appendix 5) that the basic generating function used by Debney, 

Kerr and Schild in fact satisfies equations ( 1 1 .15) and. (l 1.1 7), although 
their Killing vector is written in a rather different form. To within a 
constant multiplying factor, one can then exactly identify the variable 
£3 with Fy , and thus the choice of the form (10.26) for V leads to all 



the vacuum solutions derived in References 41 and 43 


A closer inspection of the solutions with lightlike Killing vector set 
forth in these references actually reveals that no non-trivial solutions 
exist when *=>< / 0. A recent paper by Debney ' establishes the gener- 
al result that any spacetime with a lightlike Killing vector must have 
= 0, thus there are, as we found, no solutions of (l2.22)-(l2.25). 


However, the formulation given here leads to a Killing vector only for 
the particular choice of W given by (10.26). It is not clear that there 
are no other possible choices that lead to vacuum solutions with a metric 
of the form (7.1 ) and non-zero complex expansion. For example, consider 
the form: 


W = W(xA^) 


This depends explicitly on the coordinates rather than on alone, but 
it is readily shown to satisfy (10.9) and ( 1 0. 16). It leads to equations 
similar to ( 1 1 . 1 ) , but with the 4-vector ocT replacing thus: 


OC. -fe * Ij -+ JS 

X ■+ (12.28) 


These equations require further study 
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13. The improper cane and the case of vanishing expansion . 


We now return to consider the cases where we have improper U, vanishing 

p or vanishing 9 , cases which we explicitly excluded in our development 

following equation (8.11 ). The cases of zero complex expansion have been 

considered from the point of view of gravitational radiation theory in a 
(53) 

paper by Kundt w ^ . 


Case I - 0=0. 

Referring to (8.12), 

rV 


Rotating back to the 
and noting the first 
transforms k to the 
to (8.21) : 


we see that M* will have the form: 

, / o o 0\ 

I M 0 0 

\ V 0 0 I 

original coordinate system, exactly as in 
row of the matrix R must be the vector k, 
form (l,0,0) , we have at once the result 


(13.1) 

Chapter 8, 
since R 
analogous 



(13.2) 


Since od = |S = 0 , we then have at once that 

= 0 


(13.3) 


Further, since the twist is zero, we know that the ray congruence 
must be hypersurface orthogonal, and thus that we can write in the 


form: 




Y 


(13.4) 



where (p is a scalar field and^is a scalar field (see Reference 34 # pages 
334 et seq.). 


Using (l3.3) and their definitions (7.14) and (7.16), we readily find: 


1» — -4 - — ' Xo|y 'fC 


(13.5) 


and thus, from (7.18), if we define S = G-^ 0 * we have: 


S - ~ -£<> Jy 


(13.6) 


By the use of (13.2), and changing the order of derivative operators, we 


find the important result: 

Ht;!,* 1 48 O (13.7) 

so that k^ satisfies the wave equation in the background space. 

Further, from (10.2), which is directly applicable, we have: 

= O (’3-8) 

and the second order field equations for V = 0 can therefore be written: 

- (•*.*), J* * AS,. 03-9) 

" = (S-fii)l.. + S ti (15.10) 

Expanding the left hand side of (13.10), and subtracting k^ times ( 1 3.9) 
from it, we find with the use of ( 13.7) the important result: 

S'&iu ~ Su — S|6 (13.11) 


which we can write 


= Su 


(13.12) 
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Thus, forming ^j±i± us i n S C 13-3) • we have: 

SiiU ' $lolo (13.13) 

so that S also satisfies the wave equation in the background space. 


Now, since A = 1, from ( 1 3.4) we must have: 

i (i5 - i4) 

and so 

= <t>i i/fa' (15.15) 

Since k is a unit vector, ( 1 3 • 1 5 ) at once gives: 

= 0 ( 15 . 16 ) 

Using (13.15) in ( 1 3. 1 2) gives us 

4«i / ^|o)|o 55 S|u (13.17) 

and we will therefore satisfy ( 1 3 • 1 7 ) at once if we choose: 

5 = 4 |. (,3 - ,8) 

Finally, to determine >2o , we note that ( 1 3*9 ) now takes the form 

"" ~ £ ^lojo (13.19) 


The solutions of this case are therefore characterized by: 

a) the scalar function ^ satisfies the wave equation ( 1 3 . 1 3 ) and {j) the 
eikonal-type equation ( 1 3 . 1 6 ) 

b) -£o satisfies the wave equation ( 1 3 . 1 9 ) with a source function ftlolo 

c) since (j) determines k by ( 1 3 - 1 5 ) and Ag by (l3.19)» the whole metric 
is defined by the scalar function 
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d) since the solutions are hypersurface orthogonal, solutions with 
rotation, such as the Kerr solution, cannot occur. However, 
solutions with wave properties are not excluded and are in fact 
strongly suggested by the form of the governing equations. As 
one would expect in view of the nonlinear nature of the field equations, 
solutions have a nonlinear constraint given by (13.16). 


Ca se II - p = 0 . 


This case is similar to Case I so far as the treatment of the first order 
field equations is concerned, and equations (l3.2)-(l3.7) still hold, 
except that we now have the added constraint: 


S = L = A = 0 


( 13 . 20 ) 


As before, (13.4) leads to (l3.15)» and ( 1 3 . 1 6 ) again follows from it. 
However, since we now have ( 13.20) satisfied, we have the constraint 


|v -= 0 


( 13 . 21 ) 


and (13.19) reduces to the homogeneous form: 

- o (13*22) 


Thus the relevant features of the governing equations in this case may be 
summarized as: 

a) the scalar function (j> satisfies the eikonal equation ( 1 3 . 1 5 ) 

b) the vector k^ satisfies the free space wave equation in the back- 
ground space 

c) the scalar satisfies the free space wave equation in the back- 
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ground space 

d) the vector k^ is defined by ( 1 3 • 1 5 ) 

e) the vector k^ and the scalar (|) are constrined by the relation (l3.2l). 


Case III - U is improper. 


We now explore the case where U is chosen to have negative determinant, 
which means that (8*1 1 ) can be written in the form: 


a = 


C ori & 

y — Aav0 "" Z&iQ 


(13.23) 


Thus, using (8.10), M 1 can be written as: 


/ 




\ 


0 0 0 

U l~Cc<& Swv& 

\ V S4v6 ltc<*9 j 

We find it convenient to split M' into two parts, a symmetric part from 

the lower 2x2 matrix, and a part depending only on u and v, thus: 

0 O \ 


(13.24) 


K <= J> 


O 1 -a*9 Swv9 
o Swvft / 


t 3 


(13.25) 


0 0 

V 0 0 

In rotating back to the original coordinate system, as in (8.13), we note 
that the second term of ( 1 3.25 ) is exactly as in the case of proper U, and 
thus we use that analysis to show that we must have the second term as 
k^j o k^. • For the first term of ( 1 3 • 25 ) > we note that it is symmetric, and 
thus after rotation remains symmetric. Further, it is of rank 1 and thus 


remains of rank 1 after rotation. Finally, the trace of the matrix is not 
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changed by the rotation, so the matrix after rotation has trace 2p. 
Noting that there exists a unique dyad decomposition of any symmetric 
matrix into the form: 

flj = 2 e ; w (,3 - 26) 

where Os-h, are the eigenvalues of M and e^ are the eigenvectors, we 
have in this case only one non-zero eigenvalue and thus can write the 
matrix M in the original coordinate system in the form: 

KUj (13.27) 


Use of the form (13.27) now enables us to describe the most important 

features of the solutions of the improper case without fully solving 

the field equations. For, consider the twist as defined by (9.2). 

Prom (13.27), we readily calculate that the twist is zero, and thus all 

solutions of the improper case admit a hypersurface orthogonal ray 

congruence. As noted in treating the case where 0 = 0 , this means 

that we can write k^ in the form ( 1 3. 1 5 ) where ^ is a scalar field. 

Solutions of this type, with non-zero expansion but zero twist, have 

(54) 

been studied by Robinson and Trautman . 


In the particular case where the solutions are independent of the time, 
(13.27) leads at once to: 


-iui - JijU * o 

and thus k^ is the gradient of a scalar field, 


( 13 . 28 ) 
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From (13.27) we have 

-4l\l ^ 

and then from (7.14) and (7.16) we have 


(13.29) 


(\ = *£o\l ~Al (l'3.30) 

L = + !(,-& (,3.3,) 

But since (L + A) = 2p-£o , from (13.30) and (13.31 ) we have 


L, ^ o^^)-£d 

A * 0 

^oii -Ai - 0 

so that is perpendicular to the gradient of 
From (l3. 27), we note the eigenvalue relation: 
“^vlj -€,j « Ci 


(13.32) 

(13.33) 

(13.34) 


(13.35) 


Now consider the second order field equations (13.9) and (13.10). In 
the stationary case these become: 

= 0 (,3 - 3S) 

fe.UJyv = (,3 - 57) 

Expanding the left hand side of ( 1 3. 37 ) and using ( 1 3 - 36 ) t we have 

*°c2 . 'Co 1 j Ai\j + Ao ~Ai\ jiy == -t S.f> & m 


(13.38) 
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However, ursine the fact that k^ is the gradient of y' , we note that 
-Aljlj = -fyljH , so that ( 13.38) reduces to: 

= f'ColL (15.39) 

Comparing ( 13.35) and ( 13.39). we see that we must have 

€; = -€»i; ('5.40) 

3ince M has only one non-zero eigenvalue in this case. However, this 
demands that either p = 0, or that 'i 0 \ i = 0 . 

If p = 0, (l3.27) requires that k be a constant vector, and thus ^ is 
a two-dimensional harmonic function, constant in the direction of k. 

If £o[i= 0, then &o is constant. Absorbing it into the constant term 
m of the metric (7.1 ), we see that the metric in this case must be of 
the form: 

3/W = >v> - f | V (13.41) 

where, since k^ is a unit vector, we must have ^ constrained by the 
eikonal equation: 

'hi =1 (15.42) 

We note that every case analyzed in this chapter has zero twist, thus 

only cases with proper U and non-zero complex expansion can lead to solutions 

having the rotational properties of the Kerr solutions. 
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14. Rotating solutions in the non-vacuum case . 

All the analysis and discussion to this point has dealt with the 
case where we are in vacuum and the governing field equations have 
zero matter tensor. However, it is natural to seek solutions that can 
apply within extended bodies of matter, analogous to the interior 

(2) 

Schwarzschild solution in an incompressible spherical body . The 
difficulties of solution then become formidable, mainly because much of 
the elegant structure of the field equations disappears when matter 
is present. Several approaches have been adopted in the search for 
" interior Kerr" solutions, that bear the same relation to the known 
exterior Kerr solution as does the interior Schwarzschild to the 
exterior Schwarzschild case. We will briefly discuss the method used 
in such searches, and then concentrate our attention on slow rotation, 
where as we shall see analytic solutions are attainable. In practice, 
"slow rotation" means having surface velocities small compared with the 
speed of light and this does not constitute a serious practical restric- 
ion for known stellar models. 

(55 56 ) 

To generate interior Kerr solutions, Jackson 9 J employs a method 

first given by Newman and Janis as a trick for generating the exterior 

( 57 ) 

Kerr solution from the exterior Schwarzschild solution . The essence 
of the method is the replacement of certain real coordinates by complex 
coordinates in a somewhat ad hoc manner. As Newman and Janis point out, 
it is not clear that such a procedure will even result in a solution of 
the field equations. Jackson uses a very similar method to produce 
solutions which are "complexified 11 forms of the interior Schwarzschild and 



• It is verified that the 
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the Newman, Unti and Tamburino solutions 
procedure leads to solutions that have several of the expected properties 
of interior Kerr solutions, hut unfortunately the results must be regarded 
as non-physical, since the resulting matter tensor is not positive defin- 
ite, with respect to tinelike 4-vectors* Thus there must exist observ- 
ers for whom the energy density appears negative, a most unsatisfactory 
situation. We have been unable to modify the prescription suggested by 
Newman and Janis and employed by Jackson so as to define a solution with 
the energy always positive, thus this line of attack does not appear to 
us to hold much promise for the generation of new interior solutions. 


Instead, we now limit our attention to cases where the rotation rate is 

(59) 

small. As Kerr suggested, and Boyer and Price subsequently verified, 
the product ma in the Kerr solution correspond^ to the angular momentum of 
a rotating body. Thus, let us consider the Kerr solution in the limit of 
small a, where we can neglect all powers of a higher than the first. 
Expanding the Boyer and Lindquist form of (12.18) in this way, we derive: 

!-*■">/*) - 1- . 

- ^ ( CAB 1 + JJ,*) - ctfdt 

1 f (14.1) 


This has a very suggestive form. It is simply the usual Schwarzschild 
exterior solution, with one additional cross-term added to it. We thus 
are led to investigate a metric consisting of an interior Schwarzschild 
form, modified by a cross-term of a form corresponding to that of ( 1 4 . 1 ) . 
Using the notation of Reference 8, Chapter 9, we write such a metric as: 

ds*- = 4* At 1 - df - S(eld + S^eUf) 


(14.2) 



where we assume that ”X , V and -fl are functions of r alone. The 
choice of this form, and the reasons lying behind it, deserve some disc- 
ussion. First, we expect that to first order in -TL , spherically 
symmetric configurations will remain spherical, since we certainly do 
not expect that the direction of rotation, involving the sign of ~TL , 
should be relevant to the shape of the body. Second, we know that for 

the spherically symmetric case, the most general time independent metric 

( 8 ) 

can be written in the form of (l4.2), without the term in -/L . The 
functions ^ and V are functions of r alone. 

Thus the metric chosen for the investigation of interior solutions has the 
main properties: 

1 ) When Jl = 0, the form is that of an interior Schwarzschild 
metric. Thus we can regard the added term of ( 1 4.2) as a 
perturbation to such a metric. 

2) The added term chosen to represent rotation has the same form 

as the term which represents the exterior Kerr solution to first 
order in a, added to the exterior Schwarzschild metric. 

3 ) The functions occurring in the metric, which will be determined 
in terms of the matter distribution, are functions of r alone, 
consistent with the assumption that the matter distribution rem- 
ains spherical for a treatment linear in -TL . 

Since the functions "X and V depend only on r, we see that they will 
be determined exactly as in the usual spherically symmetric Schwarzschild 
case, and cannot depend on -A. , but only on the matter distribution that 
we assume. We then expect to have an equation that JT must satisfy, in 
order to be consistent with the matter distribution and with the functions 
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'X and V • 


Suppose now that the matter is assumed to be a perfect fluid, so that the 

( 8 ) 

energy-momentum tensor has the form: ' 

~ ^ (u*Uf? ( 14 . 3 ) 

where is the 4-velocity, with index lowered using the metric of ( 1 4 . 2 ) . 
Then the 4-velocity is given by: 

= (?. jM 1 , 0, 0 , u)) (,4.4) 

where we define 10 = and assume that this is small, and comparable 

in size with J~L . The field equations that we require for this case are 
then (again see Reference 8, Chapter 9) : 


K-f = - U ( lx ) (,4 - 5) 

The evaluation of the Riemann tensor in contracted form from (l4.2) is a 
laborious but straightforward procedure, made somewhat easier by the fact 
that we can neglect quadratic or higher terms in , We find that the 
diagonal terms of the contracted Riemann tensor are exactly those of the 
usual interior Schwarzschild solution, so that we can at once assume and use 
several of the standard formulae that apply to that problem. In particular, 
we have from the diagonal terms of the field equations, the following: 


? 7 vf = ° 4-6) 

S1T(j = '/fL - ^ 4 ^ 1 (14.7) 

Sit(> = 4r°'['t r v'V - v" 

- '/*( v '- v )/+] 


(14.8) 
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The off-diagonal terms of the field equations ( 1 4 • 3 ) are identically 
satisfied to first order in jTl , except for the term R t which 
serves, with the assistance of (l4.6) - ( 1 4 *8) , to provide a single 
scalar differential equation for -fL in terms of the pressure and den- 
sity. Since all the variables are functions of r alone, this serves to 
confirm that the conjecture made earlier, that ( 1 4 • 2 ) represents a suit- 
able form for an interior metric to first order in the rotation rate. 


The governing differential equation for is: 

-e"' f- Ji" - S-rJi 1 -si * -Sl-t-v * b\i? ) (Sji 1 -Xf-SL 

L z l^. 

+ - (^4 j,) W l (14.9) 

If we now define m(r) to be the mass of the body enclosed within the 
radius r from its center, then ( 1 4 • 9) can be written in a form that 
depends only on the pressure p, density p , and mass m, the variables 
"X and v having been eliminated using the relations ( 1 4 . 6 ) - (l4.8), 



(' - 4 . - Snlf+fifjl 1 

= % n ( -ft- - u>) ( f + {>) 


The procedure is then in general the following: from the assumed equation 

of state for the fluid, the equations of Tolman, Oppenheimer and Volkoff^^ 


or some other equivalent set of relativistic equations are used to solve 
for the configuration of the body, and hence the mass distribution m(r). 
Next, the equation ( 1 4 . 1 0) is solved for the variable -fl , for an assumed 


matter rotation U) . Finally, the boundary conditions are applied at the 
edge of the rotating body, to couple to an exterior solution and hence to 
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conditions at spatial infinity. It is important to remember that the 

calculation of the mass distribution from the equation of state is done 

assuming spherical symmetry, so that a great deal of previous analysis of 

such problems is available, and provides inputs to the problem with rotation. 

When U) = constant, the entire problem can be solved in closed form for 

(29 31 ) 

the case where the equation of state is that of any perfect gas ’ . 

For the remainder of this section, we will discuss a different case, which 
corresponds to the best-known form of the interior Schwarzschild solution 
where the fluid is assumed to be incompressible . Several other cases, 
corresponding to matter shells, have been studied by Cohen and Brili^'^ 
who also point out that the function -H. is a measure of the "dragging" of 
the metric, i.e., a measure of the rotation of the inertial frame induced 
by the rotation of the body. 


For the incompressible case, we have the relations 


( 8 ). 


1 - 

where we define 

't = (l -+Vfc‘) , F = Vt* f 

= 'f, (-t c ) 


Also, solving (14.10) for LU in terms of fl , we have: 


-1 U- 


(14.11) 

(H.12) 


(14.13) 

(14.14) 


wW = H--TL' -4- -n- 


(14.15) 
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which by virtue of (l4.1l) and (14.12) can be written: 

>2 


io 


1 a + > 


(14.16) 


At the boundary of the body, where r = r^ , the metric must join in 
a continuous way to the exterior Kerr metric, expanded to first order in 
a , thus we must have from the use of ( 1 4 • 1 ) and ( 1 4. 2 ) the relation: 


_/l = - Xnms* ^ 3 


(14.17) 


(59) 

We define the total angular momentum of the body by: 


vT" ' — 'Tn./CL = y~i U 0 ) (Si. 


(14.18) 


Finally, let us consider the boundary conditions that must be applied at 
r = Tq . In order that the rotation JIl should remain finite at the 
boundary, we require that JX and Jl! should be continuous there. We 
also require, for solutions of physical interest, that P and -/X 
should be finite at r = 0 , and that w(r) should be monotonic in the 
range (0,r^) . (We do not want a star to exhibit a different sense of 
rotation at different depths within it). 


These boundary conditions, taken in conjunction with (14.16), place some 

restrictions of the form of Y1 . Thus, if jX = constant, we see at once 

from the requirement that -fll be continuous at r = r that J and hence 

0 

JX must be zero everywhere. Also, if JX contains a linear term in r, 
we find from ( 1 4 . 1 6) that u> becomes infinite at the origin. Finally, 
if we suppose that JX can be expanded as a power series in r, and that 
there is no quadratic term, use of (14.16) tells ns that we must have 
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(j(o)^/l(o), which is the case of perfect dragging, in which the inertial 
frame at the origin i3 rotating at the same speed as the star center. 

This provides an unreasonable limitation on the solutions, thus we will 
require that -fl should contain a term quadratic in r. 

kor practical application of this method, we would like (J to be not only 
monotonic in (0,r o ), but slowly varying also. Thus, for a given stellar 

density and radius, we wish to choose a functional form for SL , thus: 

oo 

-/I - + (14.19) 

-ws-a 

in which the constants c fi are chosen so as to give an acceptable slowly 
varying monotonic function 0J(r). 
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To illustrate this let us choose a fourth order function, which we 
write for convenience as: 

^ ' \ 4 | 


1 - b 


'fit ■ 


fi(r) = £l(o) 

If we match ft and ft' at the stellar boundary we obtain 


(14.20) 


b = J = ft(o)r (14.21) 

5+7x 0 V 5+7t / 

that is, the two conditions determine b and J in terms of the 
arbitrary parameter x; ft(o) plays the role of a scaling factor. The 
expression for co(r) becomes 


u)/ft(o) = 1 


9 j 

f T> 

I 2 _ 6t f-L.) 

2(5+7x) 1 

l r oV 

! 5+7x [tJ 


(14.22) 


+ ^ ~ r ) r o 5 + 14x /_r_\ 2 " 

(l + -J-) 2^ |_5+7 t 5+7t \tJ _ 

The ratio r 0 /2m^, is a convenient parameter describing the density of 
the star; for a typical neutron star it is about 5. In figure 14.1 
we have plotted u(r)/ft(o) for the cases x = 0 (no 4th order term) 
and x = .1. For x = 0,io decreases by about 20% between r = 0 and 
r = r , while for x = .1 the decrease is only about 4%, with 
u>(r) - const. =4.7 ft(o). 

The power series has been investigated also for % - 0,2,4, and 6, 
and as expected the variation in tu is extremely small. 63 



Appendix 1 . An alternative approach to the solution of the first 


order field equations . 


In Chapter 9, we mentioned an alternative method of seeking solutions of 
the first order field equations. We develop here such a method, and confine 
ourselves to the stationary case, where as we remarked in Chapter 9 the first 
order field equations serve to determine solutions completely, and where the 
analysis that follows is particularly simple and elegant. 


Prom (8.2l), in the stationary case where = 0 , we have 

"^j = — C ( Sij — •+ f £ (Al.l) 

Using this form, and calculating the quantity |J *^|> * we 

find: 

-dp ~ -fit + p ( — £i#j) (ai.2) 

This can be written concisely in terms of the complex scalar Y , thus: 


Multiplying by , we then have: 

($it — -f^) — dc — i- Eiji -$£.] Mx-j (ai.4) 

Now let us define the function F^ by: 

*” l (A1.5) 

where is a function of k^ alone. Further, let us use the equations 
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to determine as a function of the coordinates DC; . 

If we require — O , then — 0 , and thus 


' (A1 .7) 

Multiplying by (which we assume to be non-zero, consistent with the 
treatment of Chapter 8), we then have: 

^ M* ^ + ^ = 0 ( A1 « 8 ) 

and so from (A1 . 4 ) we must have: 




d-bj = 0 


(ai.9) 


In order that these differential relations in the 4 ' should be satisfied, 

<1 


we must have that: 


Y 55 . + ( ?l * _ c = 0 


(A1.10) 


which thus provides us with a relation from which we can determine Y . 
However, noting that we must also have 

- 0 (A1.1l) 

—3k 

since k is a unit vector, we see that we cannot specify three 
functions . Rather, we can choose two functions, say (j) and cj>^ » 

and then the equations that we must solve for V become: 
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^ ^ 1 ^ljb -$b 

d&i P 

<Jitj 


0 


(A1.12) 


Then k determined from (A1.6) and ^ (and thus -£•„') determined from ( A1 .12) 
completely specify the metric, as can be seen from ( 1 2 . 1 ) - (12.3). If we 
choose complicated forms for ( Pi in (A1.5), the solution for"k will be 
very cumbersome. However, for simple (|>. this method provides a simple 
alternative to the methods developed in the main body of the discussion. 

For example, let us take the simplest case, Cj). = 0 . Then from (A1.5) 
and (A1.6) vte have at once that 



= *</+ 


(A1.13) 


Also noting that we then have: 

V 

we find for the determinantal condition: 

1 -*(■&,*, '(4;z s -ix 3 )Y 

~ i-fcj 



(-€, 3(3 + ix 2 )y 

*f b 'hx 


+ i/-^3 




(A1.14) 


(A1.15) 


0 


Using (A1.13)» we see at once that both the real and imaginary parts are 
satisfied if we have Y ~ ^ t since in this case all imaginary terms 
cancel and the determinant of the remaining real coefficients is zero. 
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Comparing the results 

= a v + •, -€» = ?t(y) = '/ t Ut. i6) 

with Chapter 12, we see that we have obtained the Schwarzschild solution. 
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Appendix 2 . A direct calculation for the case where the Killing 
v ector has a vanishing time component . 

We remarked in Chanter 11 that an important relation, given by- 
equation (l1.6), wan valid even when = 0, although the result 

v;as derived explicitly assuming that ^ 0 • We demonstrate this 

assertion here. 

When = 0, from (10.23) we have 

= 0 (A 2 . 1 ) 

Also, from (8.22) we have 




(A2.2) 


We will now use the information provided by (A2„l) to calculate 
from (A2.2). The direct approach would be to use (A2.l) and two of the 
three equations of (A2.2), to give us three linear equations to be solved 
for o This becomes very messy and has to be performed for each value 

of i as a separate calculation. However, a much more elegant method is 
available, as follows: 


Since (A2.l) holds, we also have: 

~ 0 U 2 . 3 ) 

Adding this to (A2.2) gives us: 

= ^ ^ - -Mj) + f ^ 


(A2.4) 
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Now, as is readily verified by multiplication, the inverse of the matrix 

C -Aj ■+ fa fa* ) 

( 4 (4 /? 1 ~~ (a) 4j> + )/? ) (A2.5) 

where we define P = 4X u -&; . 

Multiplying (A2.4) by the inverse (A2.5), we have: 

*:i h = ( S 4f - /f) 

+ fity-U-Eiit*} 4i4 t /?) (A2 . 6) 

Now noting that = - a 1 , from ( 1 1 . 5 ) * and that P here is the same as 

that of (l0.25) when a^ = 0 , we see that (A2.6) is exactly (li.6) with 
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ArreMjjc 'i . Derivation of the gradient of the complex scalar Y . 

In Chapter 11 we wrote the formula for the gradient of without 
providing any derivation. Thin is given here. First, from (9.1 ) we 
have: 

-S ^ ~ - -fan c ( A3 . i ) 

— > <wl 

In taking the gradient of this, we will use the fact that a and k are 
independent 3-vectors, and write the general form: 

' $ -4 ^ J5 /<£ -t- C/dv-ft (A3.2) 


Multiplying by k^ , and using (9.12), which from ( 1 0.9 ) we note can be 
written as: 


C5 ^ i -& X 



(A3.3) 


we see that A can easily be found if we know B and C. Thus, in writing 

all expressions for the gradient of ^ , we will drop all terms in k . , 

3 

and recover these later by use of (A3.2). Introducing the notation " 
to mean "equal to modulo k^ ", we then have using ( 1 1 . 6 ) : 


= -flcuij = (fiiij)ii, 

— -i- Tf,iZi jeii - iftijiii 

+ Tp -itu + f Z;u 

and thus : 

?a, j - s + 7fuU it .4i tip-toi 

~ fc.ji it + f( it dSjli (».5) 
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We now use (11.6) again to evaluate the terms involving , k^^ , 

k , . , and k... , in each case dropping terms in k . . This gives us: 

r{i jji J 


^ (A3.6) 

?P Uji l£\ l — T ZfiiH *&Ac ~ T (A3.7) 

A 4 , -&ffC **“ To* p A-f- + "f|3*Ay (A3.8) 

£i<4 ^ ^^-Cr "* T ^ (A3.9) 

Using the above forms in (A3.5) then leads to: 

1«^ 4* 4 2.°^ p (A3.10) 

Now multiplying by k. and using (A3. 3) gives us the coefficient A, as: 

^ “ °^l» ““ (A3.1l) 

Thus finally from (A3.2), (A3.10) and (A3.1l) we have: 

<5y = (A3t2) 

■4 ^>/c £vje4U -4 <2. jp £/*e 


This can be written in the vector form: 

~ ( c^fo ^ +• (5 4 ~p A )-/$ 

+ ( 4 x ^3 + (iyff) 

which is exactly equation (11.8). 
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Demonstration that the complex scalar Q satisfies 
the wave equation. 


In Chapter 1 1 we stated without proof that the complex function Y 
satisfies the wave equation. To establish this result, we first write 
( 11.10) in component form: 

^ X *a + i ) + V| 




(A4.1) 


Differentiating this, we have: 

v I i. ( &.1 — -/£v /Ct + 1 £• Cj l ^ 

+ X* ( — -&IU *z° + l £-ye -&j/i ) 

"f" 0| o D Cljc 

** J 

(A4.2) 


To simplify this expression, we consider the reduction of each term. 
First, putting = P in (9.15) gives us: 



(A4.3) 


which we write: 




(A4.4) 


Also, using (11.15) and putting S='/o 


, we find: 
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Using (A4.l) to substitute for one of the Y|£ in the left hand side of 
(A4.5), and applying (A4.4), we find: 

2 Y (/Qi — -RlX^ + l £, 'P-l) 

= 3.U° 2 -ZA)V + < 27 ??,. 


hi, 


Next, using (A4.l) and differentiating to form Y|;| 0 

Yuu-k = X ? - a”) + X, 

+ i Ctjt, 4; 4j |. Zi. - i 4; 4j-|. ?, t 

Forming -*: u using (8.21) gives: 

^-j) 0 ~ 4j | ; 

+ Xlj-it 

and using this in (A4.7) then leads to: 


(A4.6) 
we find: 


(A4.7) 


(A4.8) 


y uio -Hi = -2TY5C,, +? 


|o|o 


-»• ■‘•? x Ct je , + 5pfcv J (^°-r) 

-ieijt-ftjK?,* -lAt-f (? le i e ) (A4.9) 

We can also use ( 1 1 . 6) to form , and this leads to: 

*a & ' c jU 

+ f - n l h°-T) 1 


(A4.10) 



Finally, using (11.6) again to fora £ij » multiplying 
and using ( A^I . 4 ) to reduce the result gives us: 

£ Jj'c -i^jU — T^Yio ■+ V* p ft 

+ «° 7 f> 1 + ^ ^ /*j 

ri f. now substitute the forms (A4#6), (A4#9)» (A4.10) and (A4*1l) 
which gives the result: 

$iiH ~ y, 0/0 + dc -#jn (f i %& 

~+ £('*° 2 -«‘£c) < X 3 -d?c2? V 
+ + Zfczl. 

To complete the reduction, we observe that (10.29) leads to: 

/3L % jo “ /Ct i iffi. 

and forming /C[i Xij, with the aid of (A4.1 ) then gives: 

AX l0 = ? x (&:&{**(« -T)) + fu l*' r) 

+ ie ij« Vg4t«; 

I'r.'is (4.1l) can be written as: 

y,* = -fp%, + ?‘f?‘ -+ Tp yV 

+ - u 3f?„ 


by , 

(A4.1l) 
into (A4.2), 

(A4.12) 

(A4.13) 

(A4.14) 

(A4.15) 


which means in turn that we can write the difference of equations 
and (A4.15) as: 


(A4.10) 
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% -fyc U l « x - K„ ) = (&■& - V *) 

+ r f ( - ? v) - ^ ?*( % & - /+ ?» ) 


■+• u "P V§ 0 


(A4.16) 


Now substituting this result in (A4.12), we have the final result, 
after combining terms: 






o 


(A4.17) 


Thus 


satisfies the wave equation. 
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Appendix 5 . The relation of the present treatment to alternative 
formulations. 


(43) 

Debney, Kerr and Schild have also considered metrics of the type 
(7.l) using a tetrad approach. Their results can be summarized as 
follows: 

The line element of the metric is written as: 


where the variables u,v, and 7J are defined by: 

+ Jz 1 u - 1 

>/?■? ^ ^ " 2- 1 
The function P is defined by: 


(A5.1) 


(A5.2) 


F" ' f 4 (‘J.Y + c) ~ y^) 

>, -v ( A 5.3 

— ((> Y + ^, )(w + Y* 2 } ) 

where Y is determined by solving the equation F = 0. ^ is any 

analytic function of Y, and P is defined by: 


y ' |>YV+/J.V4^.Y+C (A5.4) 

Noting that the first two terms of (A5.l) are merely a flat space metric, 
comparison with the form (7.1 ), and use of b* shows that 

for consistency of the two approaches we must require: 



(A5.5) 


We confirm this conjecture, by establishing that F^ satisfies the two 
equations ( 1 1 .15) and ( 1 1 .17). 
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First, since Y is a function of the coordinates, we note that we must have 


o = & = r Y j y -i- + — Y lj>y 

- (^y^)du '-'y^ y (A5.6) 


which we write ao 

F *1/ = -3F.^3 

where we define 

T = 


(A5.^) 


(A5.8) 


* (%Y+c,-Y(t.yt^),-((.Yf 1 ) 1 -%Yt4 

Directly from (A5.8) we have, where D denotes differentiation at constant 


Y, the relation: 

2? 5£ + if 

and thus using (A5.7) 

2& + H£y 


^ 0 


(A5.9) 




(A5.10) 


Finally, changing from the variables to z ' we find 

result: 


fDF v f - 

where (U'H is a real constant, since p and c are real, 


(A5.11) 


In an exactly similar way, we can calculate the second derivative of F. 
using the variables )W,hT . 


Y / 
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We find: 

^ + § 5 ,] = ( (A5 -’ 2) 

which in terms of the variables t,x,y,z gives: 

F, O i Fy = U5.13) 

Thus we have confirmed that Fy satisfies the same governing equations 
as U> , with the constants arising from the Killing vector written in 
rather different form. This completes the demonstration that we can 
take U> 5 F y , relating the two approaches to solving the field 
equations. 
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Fig. 1. Rotation rate as a function of radius. 


